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Abstract 



Q I Inclusive jet cross sections are measured in photoproduction at HERA using the HI detec- 

' tor. The data sample of e^p ^ e"*" + jet + X events in the kinematic range of photon vir- 

tualities < I GeV^ and photon-proton centre-of-mass energies 95 < W^p < 285 GeV 
^ I represents an integrated luminosity of 24.1 pb~^. Jets are defined using the inclusive k± 

' algorithm. Single- and multi-differential cross sections are measured as functions of jet 

^ ■ transverse energy E^'^^ and pseudorapidity i]^'^* in the domain 5 < E^'^^ < 75 GeV and 

— 1 < T]^^* < 2.5. The cross sections are found to be in good agreement with next-to- 
leading order perturbative QCD calculations corrected for fragmentation and underlying 
event effects. The cross section differential in E^^*, which varies by six orders of magni- 
tude over the measured range, is compared with similar distributions from pp colUders at 
equal and higher energies. 
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1 Introduction 



At HERA, the interaction of protons with quasi-real photons emitted from the electron^ beam 
can result in the production of jets [1, 2], for which two types of process are responsible in 
leading order (LO) quantum chromodynamics (QCD). The photon may interact as a pointlike 
particle with a parton carrying a fraction Xp of the proton momentum, in so-called direct pro- 
cesses (Fig. 1(a)). Alternatively, in resolved processes (Fig. 1(b)), the photon may develop a 
hadronic structure so that a parton carrying a fraction Xj of the photon momentum interacts with 
a parton in the proton. Due to confinement, the partons emerging from the interaction fragment 
into jets of colourless particles. The hadronic final state also includes the proton remnant and, 
in the case of resolved processes, the photon remnant and additional particles resulting from 
possible remnant-remnant interactions (the underlying event). 

The main purpose of this paper is to present inclusive jet cross sections measured using the HI 
detector and to compare them with fixed order parton level QCD predictions. After correcting 
the data and calculations to the hadron level, these comparisons offer a means of testing the va- 
lidity of perturbative QCD predictions, including the description of the partonic structure of the 
photon and the proton in terms of parton distribution functions (PDFs). The QCD-inspired mod- 
elling of non-perturbative effects in hard hadronic photoproduction can be tested as well. The 
most accurate theoretical predictions have been calculated up to next-to-leading order (NLO) in 
perturbative QCD. In order to compare these predictions with jet cross section measurements, 
the jet definition must meet certain requirements, such as infrared and coUinear safety and min- 
imal sensitivity to fragmentation and underlying event effects. The k± algorithm, originally 
proposed in [3], satisfies these requirements. 

Inclusive jet measurements at hadron colliders [4-8] and at HERA [9-13] have often been im- 
portant in the development of the understanding of QCD. In this paper, the first HI measurement 
of inclusive jet cross sections in photoproduction using the inclusive k± algorithm [14] is pre- 
sented. Compared with the last HI inclusive jet measurement [11], the integrated luminosity 
has been increased by a factor of 80 and the jet transverse energy range has been extended, now 
reaching from 5 to 75 GeV. In order to allow cross-checks with previous HERA measurements 
and comparisons with data from other colliders, the measurements have also been performed 
using a cone algorithm. 

The paper is organized as follows. In section 2, the motivations for this measurement are de- 
tailed and the phenomenology of inclusive jet photoproduction is presented. A brief description 
of the HI detector and details of the analysis procedure are given in section 3. The measure- 
ments of single- and multi-differential inclusive jet cross sections as functions of jet transverse 
energy E,p* and pseudorapidity^ rj^^^ in the laboratory frame are presented in section 4. The 
results are compared with LO and NLO QCD calculations and with inclusive jet measurements 
at pp colliders. The final section provides a summary of the results. 

'in the data taking periods used for this analysis, HERA was operated with a positron beam. However, the 
generic name "electron" will be used here to mean both electron and positron. 

^Pseudorapidity is defined asi] = — ln(tan 9/2), where 6 is the polar angle, in the coordinate system centered 
at the nominal interaction point, with respect to the positive z axis along the outgoing proton beam direction. The 
outgoing proton (photon) direction is also called forward (backward) and corresponds to 77 > (?7 < 0). 
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Figure 1 : Example LO QCD diagrams for inclusive jet photoproduction in direct (a) and re- 
solved (b) photon interactions. 

2 Inclusive jet photoproduction 

Inclusive jet cross sections are obtained by counting the number of jets found by a jet algorithm 
in a given kinematic range. The inclusive k± algorithm [14] is a modified version of the exclu- 
sive k± algorithm [15] where all hadronic final state particles are clustered iteratively^ into jets 
according to their separations in (77, 0) space"*. Here, the algorithm is applied in the laboratory 
frame. The separation parameter between jets in (i], (p) space is set to D = 1, as in [17]. An 
Et weighted recombination scheme [18], in which the reconstructed jets are massless, is used 
to maintain invariance under longitudinal boosts. To allow comparisons of the results of this 
study with previous measurements in photoproduction and with other collider data, the com- 
plete analysis has also been performed using a cone algorithm [19] with a cone radius R = 1. 
The cone algorithm has been shown to give larger hadronisation corrections than the inclusive 
algorithm (section 2.4), as is the case for jets in deep inelastic scattering (DIS) [20]. 



2.1 IMotivation 

In this paper, inclusive jet cross sections are measured over a very large E^^^ range. In the high 
E^^^ region, the high transverse momentum of the outgoing parton provides a hard scale which 
allows reliable cross section calculations to be made in perturbative QCD. It also ensures a 
reduced influence of less-well understood soft processes (fragmentation and underlying event). 
Jets at high E^^^ thus provide the most direct insight into photoproduction at the parton level. 

In the region of low E^^^, the NLO and higher order terms as well as corrections from the parton 
to the hadron level become more important, since the strong coupling as increases with decreas- 
ing scale. In the absence of a fundamental understanding of non-perturbative processes, the 
comparisons between data and theory necessarily involve phenomenological models. Matching 
the theoretical predictions with the experimental measurements at low E^^^ thus represents a 
further important test of QCD-inspired phenomenology in jet photoproduction. 

Jet photoproduction cross sections are directly sensitive to the gluon as well as the quark content 
of the photon and the proton. The proton PDFs are precisely determined [21-23] from structure 

^For more details, see e.g. [16]. 

'^(j) is the azimuthal angle in the transverse plane. 
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function measurements [24, 25]. An exception is the gluon distribution at high Xp [26]. There, 

jet photoproduction measurements [17, 27] can provide complementary information. For the 
photon, the quark density at medium and high is constrained [28-30] by F2 measurements 
at e+e" colliders [31], albeit with larger uncertainties than in the proton case. Since boson- 
gluon fusion is a higher order process compared with photon-quark scattering at e+e~ colliders, 
the gluon density in the photon is even less well constrained. Furthermore, the photon PDFs do 
not obey strict momentum sum rules, so that there is no indirect constraint on the gluon density 
in the photon. In photoproduction at HERA, higher scales can be reached than at LEP because 
of the higher centre-of-mass (cms) energy and higher 7p compared with 77 luminosity. Jet 
photoproduction cross section measurements [13, 17,27,32] thus access a largely unexplored 
domain of photon structure. 

The measurement presented here probes a wide range of E^^* and ^ quite similar to the range 
explored in a recent dijet cross section measurement [17]. Although dijet events offer better 

constraints on the hard scattering kinematics, inclusive jet measurements offer the experimental 
advantages of greater statistics and increased kinematic range and the theoretical advantage of 
naturally avoiding infrared-sensitive regions of phase space [33]. 



By considering the electron as a source of quasi-real photons of virtuality^ Q'^ and energy 
and using the QCD factorisation theorem and a factorisable jet algorithm, the electron-proton 
cross section for producing N jets {N > 2) in direct photoproduction can be expressed as: 



Here, y = E^/E^ is the fraction of the electron energy E^ carried by the photon, fy/e{y, Q"^) 
is the photon flux associated with the electron [34], fi/p{xp, /i^) is the proton PDF of parton i 
evaluated at the factorisation scale Hp and a(7 i ^ N jets) is the parton-level cross section for 
the direct subprocess ^ N jets. The cross section a is proportional to aem{lJ^R}oi^~^{lJ''R} 
at lowest order and can be expanded in powers of multiplied by perturbatively calculable 
coefficient functions, both of which depend on the renormalisation scale jiR. The kinematic 
domain over which the cross section is integrated is denoted Vt. 

Similarly, the cross section for resolved photoproduction can be written: 



Compared with Eq. 1 for direct processes, the resolved cross section in Eq. 2 includes in ad- 
dition the photon PDF of parton j, fj/j{xj, //^), evaluated at the factorisation scale //-y. Due 
to the splitting 7 qq, the QCD evolution equations of the resolved photon PDFs [35] dif- 
fer from those for the proton and lead to large quark densities at high x^, which increase with 
fij. The cross section d-{ij — > jets) describes the resolved subprocess ij^N jets and is 
proportional to ag (/i|.) at lowest order. 

^The photon virtuality is = — (fc — k')"^, where k (k') is the 4-vector of the incoming (outgoing) electron. 



2.2 QCD calculations 




(1) 
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The distinction between the direct (Eq. 1) and the resolved (Eq. 2) components is only unam- 
biguous at LO, whereas beyond LO their relative contributions depend on the factorisation scale 
m . The inclusive cross section for the production of a jet in a given kinematic range o-ep-»e+jct+x 
is obtained by summing all calculated contributions of Eqs. 1 and 2 over N, weighting by the 
corresponding number of jets inside this kinematic range. 

The partonic cross sections a in Eqs. 1 and 2 have been calculated at LO and NLO in QCD 
by several theoretical groups [33, 36-38]. These calculations differ mainly in the treatment of 
infrared and coUinear singularities. In this paper, measurements are compared with the LO and 
NLO calculations of [33], based on the subtraction method as implemented in a Monte Carlo 
program [39]. This program generates weighted parton kinematic configurations used as input 
to the inclusive k± algorithm. These calculations differ from those obtained with the phase 
space slicing method [38] by less than 2 % [40]. 

The CTEQ5M [21] parameterisation of the proton PDFs was used for the calculations. To test 
the dependence of the NLO cross sections on the choice of proton PDFs, MRST99 [22] and 
CTEQ5HJ [21] were also used, where the latter has an enhanced gluon distribution at high 
Xp. The renormalisation group equation to 2-loop accuracy was used for as and the value 
of Aqcd was chosen to match that used in the evolution of the proton PDFs (for five quark 
flavours, asiMz) = 0.118 for CTEQ5M and CTEQ5HJ, as{Mz) = 0.1175 for MRST99). 
GRV-HO [29] was chosen as the standard parameterisation of the photon PDFs. The AFG- 
HO [30] parameterisation was also used to study the dependence of the results on the choice 
of photon PDFs. The renormalisation and factorisation scales were defined as the sum of the 
transverse energies of the outgoing partons divided by 2. These scales were varied by factors 
0.5 and 2 in order to estimate the uncertainty corresponding to the missing higher-order terms. 



2.3 IMonte Carlo models 

Simulated event samples have been used to correct the data for detector effects (section 3.3) and 
to estimate hadronisation effects for the QCD calculations (section 2.4). Direct and resolved 
photoproduction events were simulated using the PYTHL\ [41], HERWIG [42], and PHO- 
JET [43] Monte Carlo generators. The generated events were passed through a GEANT [44] 
simulation of the HI detector and the same reconstruction chain as the data. 

All programs generate partonic interactions using the Bom level QCD hard scattering matrix 
elements, regulated by a minimum cut-off p on the common transverse momentum of the two 
outgoing partons. For PYTHIA and PHOJET (HERWIG), the strong coupling constant as was 
calculated by solving the 1 (2) loop renormalisation group equation using Aqcd = 200 MeV 
for 4 (5) quark flavours. GRV-LO parameterisations of the proton [23] and the photon [28] 
PDFs were used. Higher order QCD radiation effects are simulated through initial and final 
state parton showers in the leading log approximation. The fragmentation process is performed 
using the Lund string model [45] as implemented in JETSET [46] in the case of PYTHIA and 
PHOJET and using a cluster model [47] in the case of HERWIG. 

For resolved photon interactions, besides the primary parton-parton scattering, additional inter- 
actions are generated in order to simulate the effect of the underlying event. Within PYTHIA, 
these are calculated as LO QCD processes between partons from the remnants of the proton 
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and the photon. The resulting additional final state partons are required to have transverse mo- 
menta above 1.2 GcV, a value which gives an optimal description of the transverse energy flow 
outside jets for the specific photon PDFs used [11]. Soft particles accompanying the hard sub- 
process are produced in HERWIG using a soft underlying event (SUE) mechanism which is 
based on parameterisations of experimental results on soft hadron-hadron collisions. The frac- 
tion of resolved interactions which are generated with an additional SUE can be varied within 
HERWIG and has been adjusted to 35 % to match the observed level of soft activity between 
jets. PHOJET, which is based on the two-component dual parton model [48] incorporates de- 
tailed simulations of multiple soft and hard parton interactions on the basis of a unitarisation 
scheme. Due to this scheme, small variations of the lower momentum cut-off for hard parton 
interactions, set here to p^"^ — 3 GeV, do not have an influence on the results of this generator. 

2.4 Hadronisation corrections 

Since the QCD calculations refer to jets of partons, whereas the measurements refer to jets of 
hadrons, the predicted cross sections have been corrected to the hadron level using LO Monte 
Carlo programs. The hadronisation correction factors, (1 + 6hadr.), are defined as the ratio of 
the cross sections obtained with jets reconstructed from hadrons after the full event simulation 
(including parton showers, fragmentation and underlying event effects) to that from partons 
after parton showers but before fragmentation and underlying event simulation. These correc- 
tions were calculated by taking the results from two different Monte Carlo models chosen as 
described in section 3.3. The uncertainty on these corrections was taken as the quadratic sum 
of the statistical error and the systematic error given by half the difference between the results 
obtained from the two Monte Carlo models. Using the k± algorithm, the corrections were 
found to be approximately 30 % for E"^"^* < 10 GcV falling to values typically below 12 % for 
E^""^ > 20 GeV. With the cone algorithm, the corrections are around 40 % for E-^f^^ < 15 GeV 
and 20 % for E^"^* > 15 GeV. The difference between the results obtained with the two Monte 
Carlo models is typically very small and at most 10 %. 

The effects of the underlying event and of the fragmentation were also studied separately. The 
corresponding correction factors, (1 -|- Su,e.) and (1 -|- Sfrag.), were determined in the same way 
as the overall corrections factors (1 -|- Shadr.)- Here, (1 + S^.e.) is defined as the ratio of the cross 
sections obtained with jets reconstructed from hadrons with simulation of the underlying event 
to that from hadrons without simulation of the underlying event, whilst (1 + Sfrag.) is defined 
as the ratio of the cross sections obtained with jets reconstructed from hadrons to that from 
partons after parton showers, both without simulation of the underlying event. By definition, 

(1 + Shadr.) = (1 + Su.e.) ' (1 + Sfrag.)- 

Low momentum hadrons from the underlying event lead to a systematic increase of Ej,^^ and 
thus of the hadron level cross section at fixed E^jf^ The Su.e. correction is always positive 
and increases as Ejf''' decreases or t]^'^^ increases. Using the inclusive k± algorithm, for 5 < 
E^^^ < 12 GeV, Su.e. varies between - 30% at r]^^^ = -0.75 and 100% at 77^'^* = 1.25. For 
E^^* > 20 GeV, 5u.e. is always below 10 %. The effect of the underlying event is partially 
compensated by fragmentation, which has a tendency to lower the cross section. In general, 
Sfrag. is ncgativc and becomes more important as Ejf"^ decreases but also as rj^^'^ decreases, in 
contrast to 6u.e.- The Sfrag. correction is around —30 % for 5 < E^'^'' < 12 GeV and around 
-5 % for > 20 GeV. 
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3 Experimental technique 



3.1 HI detector 

A detailed description of the HI detector can be found elsewhere [49]. Here only the compo- 
nents relevant for this measurement are briefly described. 

The ep luminosity is determined with a precision of 1.5 % by comparing the event rate in the 
photon detector, located at z = —103 m, with the cross section for the QED bremsstrahlung 
process ep — > ep7. The scattered electron may be detected in the electron tagger (ETag), 
located at z — —33 m. Both detectors are TlCl/TlBr crystal Cerenkov calorimeters with an 
energy resolution of 22 

The central tracker (CT), which covers the range \r]\ < 1.5 is used to measure the trajecto- 
ries of charged particles and to reconstruct the interaction vertex. It consists of inner and 
outer cylindrical jet chambers for precise position measurement in the transverse plane, z- 
drift chambers for precise z measurement and proportional chambers for triggering purposes. 
The transverse momentum of charged particles is reconstructed from the curvature of tracks 
in the homogeneous magnetic field of 1.15 Tesla along the beam direction, with a resolution 
(t{pt)/pt ~ 0.6%pT/GeV. 

The finely segmented Liquid Argon (LAr) calorimeter [50] surrounds the tracking system and 
covers the range —1.5 < < 3.4 with full azimuthal acceptance. It consists of an electromag- 
netic section with lead absorbers and a hadronic section with steel absorbers. The total depth 
of the LAr calorimeter ranges from 4.5 to 8 hadronic interaction lengths. The energy resolution 
determined in test beam measurements is a {E)/E ^ 50 %/y/E/GeV® 2 % for charged pions. 
For jets with Ej,^^ above 20 GeV, the jet energy calibration agrees at the 2 % level with the 
Monte Carlo simulation, as determined by the Et balance in neutral current (NC) DIS and jet 
photoproduction events. At lower Ej,^*, the absolute hadronic energy scale is known to 4 %. 
The absolute resolution in rj^^^ is approximately 0.05 at = 5 GeV and better than 0.02 for 
E^,^^ > 20 GeV. 

The region —4.0 < 77 < —1.4 is covered by the SPACAL lead/scintillating-fibre calorime- 
ter [51]. Its absolute hadronic energy scale is known to 7 %. 

3.2 Event selection and reconstruction 

The data used in this paper were collected in 1996 and 1997, when electrons of energy E^, — 
27.5 GeV collided with protons of energy Ep = 820 GeV, resulting in an ep cms energy of 
300 GeV. For measurements in the region E^*^* > 21 GeV ("high" Ej,"''), the full data sample 
representing an integrated luminosity of 24.1 pb~^ was used. In addition to some activity in the 
central region, as seen by the CT, the trigger conditions required high transverse energy deposits 
in the LAr calorimeter (jet triggers). In the region 5 < i?^""* < 21 GeV ("low" E^*), where 
jet triggers suffer from proton beam-induced background, a trigger based on scattered electron 
signals in the ETag was used instead. This trigger was operated during a 'minimum bias' data 
taking period corresponding to an integrated luminosity of 0.47 pb^. The events from this 
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subsample were required to have the scattered electron detected in the fiducial volume of the 

ETag, with an energy in the range 9.6 < E[ < 19.3 GeV. The ETag geometrical acceptance, 
which is corrected for on an event-by-event basis, is always greater than 10 % in this range. 
The detection of the scattered electron ensures an improved measurement of y and hence of the 
photon-proton cms energy W.yp = AyEeEp, but reduces the available number of events by a 
factor of approximately ten, due to the restricted y range and to the limited acceptance of the 
ETag. The ETag events were also required to have no energy deposited in the photon arm of the 
luminosity system. This condition suppresses background from high rate Bethe-Heitler events 
in random coincidence with proton beam- induced background events which give activity in the 
interaction region. It also reduces QED radiative corrections. 

An interaction vertex, reconstructed from tracks in the CT and located within 30 cm of the 
nominal z position of the interaction point, was demanded. Energy deposits in the calorimeters 
and tracks in the CT were combined, in a manner that avoids double counting, in order to 

optimize the reconstruction of the hadronic final state [52], from which W^p was derived [53] 
for the "high" i?^"^* analysis. The inclusive jet sample was then defined by keeping all events 
for which at least one jet was reconstructed with the inclusive k± algorithm in the kinematic 
domain: 

-1 < V^^* < 2.5 ; E^^^ > 21 GeV ("high" E^*) ; (3) 

-1 < rj^^* < 2.5 ; E^^^ > 5 GeV ("low" Ef'). (4) 

The rj^'^^ range was chosen to ensure that the jets were well contained in the LAr calorimeter. 
For "high" Ej,^* events, the kinematic region was restricted to 

g2 < 1 GeV^ ; 95 < W^p < 285 GeV. (5) 

The Q'^ range is implied by the absence of the scattered electron in the main HI detector. For 
"low" E^^* events, the tagging of the electron already restricted the kinematic region to 

Q2 < 10-2 GeV^ ; 164 < W^p < 242 GeV. (6) 

A number of requirements were made in order to suppress the non-ep background. For "high" 
E^j.^'^, the vertex condition was sufficient to reduce the contamination from beam-induced back- 
ground to a negligible level. Background originating from cosmic showers and beam halo 
muons was rejected using a set of topological muon finders [54]. In addition, the total missing 
transverse momentum ^j, was required to be small compared with the total transverse energy 
Et by applying the cut ^j,/^/E^ < 2.5 GeV^. The overall non-ep background contamination 
in the "high" Ej,^^ sample was then estimated to be below 1 %. For "low" E,p* events, since 
W^p can be measured in two independent ways (using either the energy of the scattered electron 
or the hadronic final state), consistency between the results of the two methods was required. By 
fitting the distribution of the z position of the vertex with the sum of a Gaussian and a constant, 
the "low" E^^^ sample was estimated to have a non-ep background contamination of about 2 %. 

Further cuts were applied to reduce backgrounds from other ep collision processes. For the 
"high" E^^* jet sample, the only significant one is the NC DIS process^, where either the scat- 
tered electron or the current jet enters the inclusive jet selection. Events with a scattered electron 

*The charged current DIS background was already completely suppressed by the ^j, cut. 
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candidate found using electron identification algorithms [25] were thus rejected. The remaining 

ep background in the "high" E"^"^* sample was estimated to be below 1 % from a study using 
a simulated sample of NC DIS events from the DJANGO [55] Monte Carlo generator. In the 
"low" Ejf"*^ sample it is completely negligible. 

The event samples finally consist of 15 388 jets reconstructed in 11 801 events for "high" E^'^^ 
and 26 848 jets reconstructed in 21001 events for "low" For both samples, the total 

background of around 2 % was subtracted. The inefficiency due to selection cuts is below 3 %. 

3.3 Cross section determination 

To obtain the inclusive jet cross section, each of the N jets reconstructed in a given range is 
assigned a weight calculated as the inverse of the event-by-event trigger efficiency e. The cross 
section obtained at the detector level is then corrected by a factor C for inefficiencies due to 
the selection procedure and for migrations caused by the detector response. For instance, the 
double-differential cross section averaged over a range AEj^^^Arj^'^^, is defined as: 

where C is the integrated luminosity. 

The trigger efficiency e was determined from data by using events triggered independently. For 
the "high" E^^^ analysis, e was parameterised as a function of the Et and i] of the leading jet and 
was always greater than 80 %, reaching 100 % for Et > 35 GeV. For the "low" E^^*^ analysis, 
e was found to depend on the multiplicity of CT tracks in the event, with a mean value of 97 %. 

Two Monte Carlo programs (section 2.3) were used to correct the data from each event sample 
for detector effects. For the "high" E^''^ sample, HERWIG and PYTHIA were chosen. A 
reasonable description of the observed energy flow around the jet axis was obtained with both 
programs, provided the underlying event or multiple interaction mechanisms were included 
in the Monte Carlo simulations [40]. For the "low" Ejf''^ sample, PHOJET, which has been 
shown to give the best description of energy flow distributions [11] and jet profiles [40], was 
chosen together with PYTHIA. The mean correction factors C calculated for each measurement 
interval with the two Monte Carlo models were found to lie between 0.9 and 1.6 for the "high" 
E^^* analysis and between 0.5 and 1.6 for the "low" Ej,^* analysis. 

The bin-to-bin migrations are important due to the steeply-falling shape of the spectrum. 
The bin widths were chosen to measure cross sections in as many intervals as possible whilst 
ensuring that stability and purity^ were greater than 30 %. These criteria restrict the r]^"^ cross 
section measurements to different r]^"^ ranges depending on the E-^"^* range considered. At 
"high" E^^*, the problematic region is that of negative 77^^*. Due to the Lorentz boost between 
the 7p cms and the laboratory frame, r]^'^^ = corresponds to a cms pseudorapidity^ rj* of about 
—2, which is well into the photon hemisphere. The cross section thus falls most steeply with 
increasing in this region. At "low" E^^^, the influence of the proton remnant compromises 
measurements towards higher 77 -^^^ 

^The stability S (purity V) is defined as the number of jets which are both generated and reconstructed in 
an analysis bin, divided by the total number of jets that are generated (reconstructed) in that bin. By definition, 

C = S/V. 
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3.4 Systematic uncertainties 



The following systematic uncertainties have been considered : 

• The uncertainty in the absolute hadronic energy scale of the LAr calorimeter (2 % for 
"high" Ep^ and 4 % for "low" E^^') leads to an uncertainty of typically 10 % for "high" 
E^^^ and between 10 % and 20 % for "low" E^^\ 

• The 7 % uncertainty in the hadronic SPACAL energy scale results in an uncertainty of 
1 % at "high" E^^^ and is negligible at "low" Ef\ 

• The 3 % uncertainty in the fraction of the energy of the reconstructed hadronic final state 
carried by tracks leads to an uncertainty of less than 1 % for "high" E^"^* and of 2 % to 
4 % for "low" E^^\ 

• The background subtraction leads to an uncertainty of 1 %. 

• The statistical uncertainty in the trigger efficiency determination leads to an uncertainty 
of 1 % or less. 

• The uncertainty in the integrated luminosity results in an overall normalisation error of 
1.5%. 

• The uncertainty in the correction for detector effects was taken to be half the differ- 
ence between the correction factors calculated from the two Monte Carlo programs. It 
is smaller than 8 % for "high" Ef^ and smaller than 10 % for "low" Ef\ 

All systematic uncertainties are added in quadrature. The resulting uncertainty ranges from 
10 % to 20 % for "high" E^'^ and from 15 % to 30 % for "low" E^^^ and is dominated by the 
normalisation uncertainty due to the LAr calorimeter energy scale. 



4 Results 

In this section, inclusive jet cross sections are presented, corrected for detector effects and mea- 
sured in different kinematic regions as functions of E^'^'^ and rj^'^* in the laboratory frame. Good 
agreement with previous data [9-12] has been found when using the cone algorithm [40]. The 
results are also consistent with recently published ZEUS data [13] using the k± algorithm. The 
numerical results using the k± algorithm (Tabs. 1 to 5) are given as differential cross sections 
averaged over the quoted ranges. Those obtained with the cone algorithm (Tabs. 6 and 7) are 
given at the average value in each analysis interval, determined according to the Monte Carlo 
simulation. The results are shown in Figs. 2 to 9. In the E^'^'^ spectra (upper part of Figs. 2 
to 4), all systematic uncertainties are added in quadrature with the statistical uncertainty and 
are shown as error bars. The inner error bars denote the statistical and the outer the total uncer- 
tainty. In all other figures, the LAr calorimeter energy scale uncertainty is not included in the 
error bars, but is shown separately as a hatched band. The results are compared with LO and 
NLO QCD predictions first, then with similar measurements at pp colliders. Unless otherwise 



9 



stated, the QCD prediction is calculated at NLO with the standard settings described in sec- 
tion 2.2 and corrected for hadronisation, as explained in section 2.4. Its uncertainty is shown as 
a shaded band divided into two parts. The inside (light) part is the uncertainty associated with 
the hadronisation corrections and the outside (dark) part is the uncertainty associated with the 
choice of the renormalisation and factorisation scales. These uncertainties are added linearly. 
When presented (Figs. 2 to 5), the relative differences are always defined with respect to this 
standard NLO QCD prediction. The total hadronisation correction factors (1 + 5hadr.) and their 
errors are given in the tables, together with the correction factors associated with the effects of 
fragmentation, (1 + 6 frag.) and of the underlying event, (1 + Su.e.)- 



4.1 cross sections 

The measured differential e+p cross section da/d£'|.'^* for inclusive jet production integrated 
over — 1 < r]^'^* < 2.5 in the kinematic region defined by < 1 GeV^ and 95 < Wyp < 
285 GeV is shown for E^""* > 21 GeV in Fig. 2 (top). The LO QCD calculation fails to repro- 
duce the normalisation of the distribution. Both NLO predictions, with and without hadronisa- 
tion corrections, reproduce the measured distribution very well. As shown in Fig. 2 (bottom), 
the uncertainty due to the renormalisation and factorisation scales is of the order of 10 %. The 
calculated cross sections using the GRV photon PDFs are typically 5 % to 10 % larger than 
those obtained with AFG. To show the sensitivity to the proton PDFs, the predictions using 
GRV for the photon and MRST99 or CTEQ5HJ for the proton are also shown. Compared with 
CTEQ5M, MRST99 and CTEQ5HJ give almost the same prediction at relatively small E^''^ 
but show differences as i?^"^* increases. The prediction using MRST99 decreases relative to that 
using CTEQ5M by 5 % over the measured Ejf''' range while that using CTEQ5HJ increases by 
8 %. Within the errors, the NLO QCD calculations with each of the PDFs choices describe the 
magnitude and the shape of the measured inclusive Ej,^* spectrum very well, up to the highest 
measured Ej,^^ values. 

In Fig. 3 (top), da/dE^^* is presented for two separate W^p intervals 95 < W^p < 212 GeV 
and 212 < W^p < 285 GcV. At higher W^p, the Ej,^* spectrum is harder and extends to higher 
i?^"^* values, as expected. The predictions of the two combinations of photon and proton PDFs 
which give the lowest (AFG for the photon and MRST99 for the proton) and the highest (GRV 
for the photon and CTEQ5HJ for the proton) cross sections are also shown. As can be seen in 
Fig. 3 (bottom), all these NLO QCD calculations describe the magnitude and the shape of the 
inclusive Ej,^* spectra measured in the two Wjp ranges. 

A measurement of the inclusive jet cross section over the whole E,p* range was performed by 

combining the "low" and "high" Erp^ data samples. In order to do this, the same Wjp cut was 
applied to the "high" as to the "low" E^""^ data sample, i.e. 164 < W^p < 242 GeV. The 
"low" E^'^'^ cross section was also corrected to correspond to the same range^ as the "high" 
E:p^ sample. The measured cross section is shown in Fig. 4. The analyses of the "low" and 
"high" Ej,^* samples agree well in their domain of overlap. The measured da/dE,p^ cross 

'The "low" E^""* cross section was multipUed by the ratio Rf = F{1 GeV^)/i^(0.01 GcV^), where F(Q^„^) 
is the integral of the photon flux ,f-y/e{y,Q^) (see Eqs. 1 and 2) over Q"^ < Q^ax the range 0.3 < y < 
0.65, which corresponds to the chosen W^p range. The numerical integration yields -F(l GeV^) = 0.0181, 
F(0.01 GeV^) = 0.0127 and Rp = 1.43. 
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section falls by more than 6 orders of magnitude between E^^ — 5 and 75 GeV and is well re- 
produced by the theoretical prediction. The NLO contribution and the hadronisation corrections 
are both needed to give a good agreement of the calculation with the measured data. 

Following a procedure previously applied to the inclusive charged particle photoproduction 
cross section [56], the power-law [57] function oc (1 + E^^^ / Eta)'^ was fitted to the inclusive 
jet cross section 1/E^^^ ■ da/dE^^* obtained from the data of Fig. 4. Good fits with stable 
results could only be obtained in the region 5 < E.p'' < 35 GeV. The fit gave the results 
Et,o = 2.4 ± 0.6 (stat.) t^l (syst.) GeV and n = 7.5 ± 0.3 (stat.) IqIs (syst.). The value of 
the exponent is in agreement with that obtained in [56]: n — 7.03 ± 0.07 (stat.) ± 0.20 (syst.). 

4.2 77^^* cross sections 

The measured differential e~^p cross section da/di]^'^^ in the range —1 < r]^^^ < 2.5 for in- 
clusive jet production integrated over 21 < i?^"* < 75 GcV, < 1 GcV^ and 95 < W^p < 
285 GeV is shown in Fig. 5. The calculated cross sections using different combinations of the 
photon and proton PDFs give a good description of the data, within the experimental and theo- 
retical errors. The normalisation of the data can only be described when the NLO contribution 
is included. The description is not significantly improved by the hadronisation corrections. 

In Fig. 6, da/dri^^* is presented in three different intervals of Ej,^^ and compared with NLO 
QCD predictions. The hadronisation corrections correspond to an increase (decrease) of the 

pure partonic prediction in the forward (backward) region. The increase in the forward region 
is due to the influence of the proton remnant leading to significant underlying energy, while 
the decrease in the backward region reflects the escape of partonic energy from the jet due to 
fragmentation (section 2.4). Within the errors, the data are well described by the NLO QCD 
predictions. 

In order to study the cross section more differentially, measurements of da/drj^^* in two regions 
of W^p and three intervals of E^^* are presented in Fig. 7. The maximum of the cross section is 
shifted towards low r/^*^* values at higher W^p due to the decreased forward boost of the hadronic 
cms relative to the laboratory frame. The NLO QCD predictions with and without (not shown 
in Fig. 7) hadronisation corrections are in general in good agreement with the measured cross 
sections. NLO QCD calculations using different combinations of photon and proton PDFs also 
give good agreement with the data. The precision of the measurement as well as the theoretical 
uncertainties do not allow any firm conclusion to be drawn on which combination of PDFs is 
favoured by the data. 

The da/drj^'^^ measurement for the "low" E^^* sample is presented in Fig. 8 for two intervals 
of in the kinematic region < 10^^ Q^y^ and 164 < W^p < 242 GeV. The NLO 
predictions are in agreement with the data in the range 12 < E^^* < 21 GeV, provided the 
hadronisation corrections are included. In the lowest E^'^^ range 5 < Ef' < 12 GeV, however, 
the agreement with the NLO predictions including hadronisation corrections is marginal and the 
data seem to indicate a rise of the cross section with increasing i]^^^ which is faster than in the 
theoretical predictions. This may be the result of a failure of the LO Monte Carlo to accurately 
estimate the hadronisation corrections. Inadequacy of the photon PDFs in this kinematic range 
or the absence of higher order corrections beyond NLO may also be responsible. 
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4.3 Comparison with pp collider results 



It is interesting to compare the present photoproduction measurements with similar pp data in 
order to see the effects arising from the different structure of the photon and the proton. The 
differential e+p cross section da/dE^'^^ was measured in the range —1 < ry-'^* < 2.5, as in 
Fig. 4, but using the cone algorithm with cone radius i? = 1 to match the procedure used 
for the only available pp data [4] at comparable cms energy ^/s = 200 GeV. The results are 
presented in Tab. 6. Up to a normalisation factor, the Ejf^'' dependence of the data is compatible 
in the region of overlap 5 < Erp* < 22 GcV with that of [4]. Monte Carlo studies show 
that the difference in cms jet pseudorapidity ranges between the 7p (—3.0 < rj* < 0.5) and pp 
(1 77* I < 1.5) data does not affect this conclusion. 

To allow comparisons with pp measurements at different energies, the scaled 'jp cross section 



where [E^^^^p^^*") is the four-vector of the jet, was measured as a function of the dimensionless 
variable xt — 2E^'^^/W^p. In the naive parton model, S{xt) is independent of cms energy for 
the same colliding particles. The differential e'^p cross section da/dE^^* was measured with 
the cone algorithm in the restricted range 1.5 < 7]^^* < 2.5 and E"^^* > 8 GeV. It was then 
transformed into S{xt) at a fixed W-yp = 200 GeV averaged over the range |?7*| < 0.5 using the 
Monte Carlo models to evaluate the correction factors and their uncertainties. 

The S{xt) distribution is presented in Tab. 7 and Fig. 9. In the figure, it is compared with data 
from pp scattering obtained by the UAl [4,5] and DO [6,7] collaborations at various energies 
using the cone algorithm. The pp data were transformed into S{xt) using the E^'''^ value at 
the centre of each measurement interval and were scaled by factors of 0{aem/c(s) such that 
S{xt) approximately matches that from the photoproduction data at xj- ~ 0.1. Despite the 
differences in the rf ranges of measurement and in the details of the analysis procedure'", all pp 
data are in approximate agreement after the scaling factors are applied. Within the experimental 
uncertainties, the shape of the S{xt) distribution for 7p is compatible with those from pp data 
in the region xt 0.2, where the resolved photon leads to a similar behaviour of the scaled 
cross section to that for a hadron except for the overall normalisation. At larger xt, the shape 
of the 7p cross section begins to deviate from that for pp. As can be inferred from the PYTHIA 
predictions for the full -jp cross section and for the contribution from resolved photons, this is 
due to the enhancement of the resolved photon quark density relative to that of the proton at 
large momentum fractions, as well as the increasing proportion of direct photon interactions. 
The direct photon contribution involves the convolution of only one set of PDFs and dominates 
the scaled cross section at the largest xt- 

'"UAl measured cross sections in the range \ri*\ < 1.5 for y/s = 200 GcV and \ri*\ < 0.7 for y/s — 630 GeV, 
using a cone radius R = 1 and no jet pedestal energy subtraction. DO measured cross sections in the range 
|?7*| < 0.5 using a cone radius i? = 0.7 and jet pedestal energy subtraction. 
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5 Summary 



A new measurement of inclusive jet production cross sections in quasi-real photoproduction 
(Q^ < 1 GcV^) has been presented, based on an integrated luminosity of 24.1 pb^^ of e+p data 
collected by the HI experiment in the years 1996 and 1997. Compared with the last published 
HI result [11] on this topic, this measurement represents an increase in luminosity by a factor 
of 80. The jets were selected using the inclusive k± algorithm in the pseudorapidity range 
— 1 < rj^^* < 2.5 in the laboratory frame. The photon-proton centre-of-mass energy range of 
the measurement for jets with transverse energies > 21 GeV is 95 < Wjp < 285 GeV. 
The measurement could be extended down to E;^'' > 5 GcV by using a sample of data with 
integrated luminosity 0.47 pb \ collected in a data taking period with a dedicated trigger. 
There, the kinematic range of measurement was < 0.01 GeV^ and 164 < Wjp < 242 GeV. 

The measured cross sections were corrected to the hadron level and compared with leading 
order (LO) and next-to-leading order (NLO) QCD calculations, with and without fragmenta- 
tion and underlying event corrections. The LO QCD calculations are unable to reproduce the 
normalisation of the experimental data. The NLO QCD calculations, using various available 
photon and proton PDFs, describe the measured distributions both in normalisation and shape 
over the whole E'^"* and r/-''^* range within the experimental and theoretical uncertainties. For 
E^'^'^ > 21 GeV, the hadronisation corrections to the NLO QCD calculations only slightly im- 
prove the agreement with the data, whereas for 5 < E^^* < 21 GeV, good agreement can only 
be obtained with the hadronisation corrections. The current precision of the experimental re- 
sults as well as of the theoretical predictions does not allow one to discriminate between the 
different photon and proton PDFs using these data alone. However, the information obtained 
from these measurements could be used to constrain the photon and proton PDFs in global fits 
of experimental results. 

To compare with previous measurements at HERA and at pp colliders, the inclusive E^"'' dif- 
ferential cross section was also measured for jets defined using the cone algorithm with R—1. 

The shape of the scaled cross section S{xt), as a function of the dimensionless variable 
xt = 2E^^^ /W^p at W^p = 200 GcV for \'q*\ < 0.5, is compatible with that of similar 
measurements for xt < 0.2. The shapes for and pp are different at larger xt, where resolved 
photon structure at large Xy and direct photon interactions become important. 
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7-1 jet 

Ej. range 


da/dEir' 




(1 + 5 frag.) (1 + K.e) (1 + ^hadr.) 


[GeV] 




[pb/GeV] 





95 < < 285 GeV 



21. 


.28 


65.4 


±0.6 


±1.6 


+5.5 

-5.3 


0.95 ±0.02 


1 


08 ±0.01 


1 


03 ±0.03 


28. 


.35 


14.0 


±0.3 


±0.3 


+1.4 
-1.2 


0.95 ±0.02 


1 


06 ±0.02 


1 


00 ±0.04 


35. 


.42 


3.56 


±0.14 


±0.09 


+0.39 
-0.33 


0.95 ±0.01 


1 


04 ±0.01 


1 


00 ±0.02 


42. 


.52 


0.908 


±0.060 


±0.018 


+0.107 
-0.095 


0.95 ±0.01 


1 


01 ±0.02 





96 ±0.01 


52. 


.62 


0.192 


±0.028 


±0.012 


+0.023 
-0.021 


0.97 ±0.03 


1 


01 ±0.03 





98 ±0.04 


62. 


.75 


0.0483 


±0.0121 


±0.0018 


+0.0063 
-0.0071 


0.96 ±0.04 


1 


00 ±0.05 





96 ±0.04 



95 < W^^ < 212 GeV 



21. 


.28 


32.7 


±0.4 


±0.7 


+2.6 
-3.0 





94 ±0.02 


1 


07 ±0.01 


1 


00 


±0.03 


28. 


.35 


6.21 


±0.19 


±0.18 


+0.60 
-0.50 





94 ±0.01 


1 


04 ±0.01 





98 


±0.02 


35. 


.42 


1.51 


±0.09 


±0.03 


+0.15 

-0.14 





94 ±0.02 


1 


04 ±0.01 





98 


±0.01 


42. 


.52 


0.236 


±0.030 


±0.008 


+0.025 
-0.024 





93 ±0.02 


1 


00 ±0.03 





93 


±0.02 


52. 


.62 


0.0360 


±0.0115 


±0.0009 


+0.0052 
-0.0041 





92 ±0.04 


1 


00 ±0.05 





92 


±0.04 


62. 


.75 


0.00511 


±0.00365 


±0.00019 


+0.00041 
-0.00115 





90 ±0.11 





98 ±0.12 





88 


±0.10 



212 < < 285 GeV 



21. 


.28 


32.8 


±0.4 


±1.0 


+2.8 
-2.3 





97 ±0.01 


1 


10 ±0.02 


1 


06 ±0.03 


28. 


.35 


7.81 


±0.21 


±0.21 


+0.83 
-0.71 





95 ±0.02 


1 


07 ±0.03 


1 


03 ±0.06 


35. 


.42 


2.05 


±0.11 


±0.07 


+0.23 
-0.20 





96 ±0.01 


1 


05 ±0.03 


1 


01 ±0.03 


42. 


.52 


0.676 


±0.053 


±0.017 


+0.083 
-0.071 





97 ±0.02 


1 


01 ±0.02 





98 ±0.02 


52. 


.62 


0.157 


±0.026 


±0.012 


+0.018 
-0.017 





99 ±0.03 


1 


01 ±0.04 


1 


00 ±0.05 


62. 


.75 


0.0434 


±0.0117 


±0.0018 


+0.0061 
-0.0057 





97 ±0.04 


1 


00 ±0.05 





97 ±0.05 



Table 1: Measured differential e'^p cross section da/dEjf for inclusive jet photoproduction 
(Q^ < 1 GeV^j, integrated over —1 < 77^*^* < 2.5 in three regions ofW^p. Jets are defined 
using the inclusive k± algorithm. The statistical (A.stat.), systematic (A.syst.) (excluding LAr 
energy scale) and LAr energy scale (A^.s.) uncertainties are shown separately. The correction 
factors applied to the NLO QCD predictions are also shown separately as (1 + 5 frag.), for 
fragmentation, (1 + 5u.e.), for the underlying event, and the product (1 + Shadr.) for the total 
hadronisation correction. 
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range 


do/dE^^ Agtat. ^syst. ^e.s. 


( 1 + 8 frag. ) ( 1 + (^«.e. ) ( 1 + hadr. ) 


[GeV] 


[pb/GeV] 





161 < < 212 GeV ; Q'^ < 0.01 GeV^ 



5. 


..8 


24600 


±200 


±1600 


+3000 
-2900 


0.72 ±0.06 


1.77 ±0.23 


1.25 ±0.06 


8.. 


. 12 


3070 


±60 


±230 


+520 
-470 


0.80 ±0.08 


1.66 ±0.11 


1.31 ±0.06 


12. 


..16 


505 


±26 


±30 


+94 
-84 


0.87 ±0.08 


1.43 ±0.10 


1.23 ±0.04 


16. 


. .21 


126 


±11 


±6 


+27 
-19 


0.83 ±0.11 


1.26 ±0.09 


1.04 ±0.09 


21. 


..28 


28.3 


±6.1 


±4.6 


+10.1 
-6.9 


0.86 ±0.17 


1.23 ±0.13 


1.05 ±0.17 



164 < W^p < 242 GeV ; < I GeV^ 



21. 


.28 


30.1 


±0.4 


±0.8 


+2.4 
-2.5 


0.95 ±0.02 


1.09 ±0.01 


1.04 ±0.03 


28. 


.35 


6.74 


±0.19 


±0.18 


+0.63 
-0.58 


0.95 ±0.01 


1.07 ±0.02 


1.01 ±0.04 


35. 


.42 


1.66 


±0.10 


±0.04 


+0.18 
-0.14 


0.96 ±0.01 


1.04 ±0.02 


1.00 ±0.02 


42. 


.52 


0.417 


±0.041 


±0.013 


+0.043 
-0.040 


0.96 ±0.02 


1.01 ±0.03 


0.97 ±0.02 


52. 


.62 


0.0773 


±0.0174 


±0.0066 


+0.0113 
-0.0119 


0.95 ±0.03 


1.00 ±0.04 


0.95 ±0.03 


62. 


.75 


0.0132 


±0.0059 


±0.0014 


+0.0016 
-0.0019 


0.94 ±0.07 


0.95 ±0.08 


0.89 ±0.09 



Table 2: Measured differential e^p cross section da/dEj.^ for inclusive jet photoproduction, 
integrated over —1 < r]^^^ < 2.5 in the kinematic region 164 < W^p < 242 GeV (see Tab. 1 
caption for further details). 



18 



fjjet range 


da / dri^''^ Astat. A^yst. A^.s. 


( 1 + Sfrag. ) ( 1 + ^u.e. ) ( 1 + Shadr. ) 




[Pb] 





21 < E^^^ < 75 GeV 



-1. 


..0 


37.8 


±1.2 


±2.1 


+3.0 
-2.5 


0.84 ±0.03 


1.04 ±0.03 


0.87 ±0.05 


0.. 


.0.5 


173 


±4 


±4 


+12 
-12 


0.90 ±0.03 


1.06 ±0.02 


0.96 ±0.05 


0.5 


..1 


257 


±5 


±7 


+21 
-23 


0.94 ±0.02 


1.06 ±0.02 


1.00 ±0.04 


1.. 


.1.5 


253 


±4 


±6 


+25 
-20 


0.96 ±0.01 


1.07 ±0.02 


1.02 ±0.03 


1.5 


..2 


237 


±4 


±7 


+22 
-19 


0.99 ±0.01 


1.09 ±0.01 


1.08 ±0.01 


2.. 


.2.5 


186 


±4 


±4 


+17 
-19 


1.01 ±0.01 


1.11 ±0.02 


1.12 ±0.01 



21 < Ef' < 35 GeV 



-1. 


..0 


37.6 


±1.1 


±2.1 


+2.9 
-2.5 


0.84 ±0.03 


1.04 ±0.03 


0.87 ±0.05 


0.. 


.0.5 


166 


±4 


±4 


+11 
-11 


0.90 ±0.03 


1.06 ±0.02 


0.96 ±0.05 


0.5 


..1 


241 


±4 


±6 


+19 

-21 


0.94 ±0.02 


1.07 ±0.02 


1.00 ±0.04 


1.. 


.1.5 


233 


±4 


±6 


+23 
-18 


0.96 ±0.02 


1.07 ±0.02 


1.02 ±0.03 


1.5 


..2 


220 


±4 


±7 


+20 
-18 


0.99 ±0.01 


1.10 ±0.01 


1.09 ±0.01 


2.. 


.2.5 


174 


±4 


±4 


+16 
-18 


1.01 ±0.01 


1.11 ±0.02 


1.13 ±0.02 



35 < E^^^ < 52 GeV 



0...0.5 
0.5... 1 
1...1.5 
1.5. . .2 
2. . .2.5 


7.56 ±0.75 ±0.64 to^l 
14.8 ±1.1 ±0.4 t\i 
18.3 ±1.2 ±0.4 tli 
15.1 ±1.1 ±0.3 til 
11.5 ±1.0 ±0.2 tli 


0.87 ±0.03 1.01 ±0.02 0.87 ±0.03 
0.94 ±0.02 1.03 ±0.02 0.96 ±0.03 
0.96 ±0.01 1.03 ±0.02 0.99 ±0.02 
0.98 ±0.01 1.04 ±0.02 1.02 ±0.01 
0.99 ±0.02 1.06 ±0.02 1.05 ±0.02 


52 < E^,^* < 75 GeV 


0.5... 1 
1...1.5 
1.5... 2 
2... 2.5 


1.16 ±0.33 ±0.08 tofi 
1.69 ±0.37 ±0.10 1^:23 
1.84 ±0.39 ±0.12 toil 
0.458 ±0.189 ±0.040 toml 


0.93 ±0.05 1.03 ±0.05 0.96 ±0.05 
0.97 ±0.04 0.97 ±0.05 0.94 ±0.06 
0.99 ±0.04 1.03 ±0.04 1.02 ±0.04 
0.98 ±0.06 1.04 ±0.07 1.02 ±0.08 



Table 3: Measured differential e~^p cross section da / drj^'^^ for inclusive jet photoproduction 
(Q^ < 1 GeV^ j, integrated over four E^^^ ranges in the kinematic region 95 < W^p < 285 GeV 
(see Tab. 1 caption for further details). 
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fjjet range 




^stat. 

[Pb] 






(l + 6f ) 


(1 + S ) 


(1 + Sh ^ ) 

\-'- ~ "naar.J 








21 < E^^^ < 35 GeV 


; 95<w^p 


< 212 GeV 




0.. 


.0.5 


32.6 


±1.5 


±0.9 


+1.7 
-3.2 


0.79 ±0.04 


1.02 ±0.04 


0.81 ±0.06 


0.5 


. . . 1 


114 


±3 


±4 


+10 
-9 


0.90 ±0.02 


1.06 ±0.02 


0.95 ±0.03 


1 . . 


. 1.5 


141 


±3 


±3 


+12 

-12 


0.93 ±0.02 


1.05 ±0.02 


0.98 ±0.04 


1.5 


. . . 2 


142 


±3 


±5 


+13 
-11 


0.98 ±0.01 


1.06 ±0.02 


1.04 ±0.02 


2 


. 2.5 


111 


±3 


±.3 


+9 
-il 


0.99 ±0.02 


1.10 ±0.03 


1.09 ±0.02 








21 <E 


< 35 GeV 


; 212 < W^p 


< 285 GeV 




-1 


...0 


37.3 


±1.2 


±1.9 


+3.1 
-2.9 


0.85 ±0.03 


1.04 ±0.03 


0.88 ±0.05 


0.. 


.0.5 


133 


±3 


±4 


+9 
-8 


0.94 ±0.03 


1.07 ±0.02 


1.01 ±0.05 


0.5 


. . . 1 


127 


±3 


±6 


+9 

-12 


0.97 ±0.02 


1.08 ±0.04 


1.05 ±0.05 


1 . . 


. 1.5 


91.7 


±2.6 


±3.7 


+11.6 
-6.2 


1.01 ±0.02 


1.10 ±0.03 


1.11 ±0.03 


1.5 


. . . 2 


78.7 


±2.6 


±2.2 


+6.8 
-6.3 


1.02 ±0.02 


1.16 ±0.03 


1.19 ±0.03 


2 . . 


. 2.5 


59.6 


±2.3 


±2.4 


+7.4 
-3.9 


1.05 ±0.03 


1.14 ±0.03 


1.20 ±0.03 








35 < Ef* < 52 GeV 


; 95 < W^p 


< 212 GeV 




0.5 


. . . 1 


2.72 


±0.46 


±0.17 


+0.16 
-0.18 


0.82 ±0.04 


1.03 ±0.04 


0.85 ±0.05 


1 . . 


. 1.5 


8.07 


±0.81 


±0.19 


+0.76 
-0.72 


0.93 ±0.02 


1.02 ±0.02 


0.95 ±0.02 


1.5 


. . . 2 


8.27 


±0.81 


±0.20 


+0.87 
-0.75 


0.96 ±0.02 


1.02 ±0.02 


0.98 ±0.02 


2 . . 


.2.5 


6.57 


±0.71 


±0.19 


+0.84 
-0.71 


0.97 ±0.03 


1.05 ±0.03 


1.02 ±0.02 








35 <E 


< 52 GeV 


; 212 < W^p 


< 285 GeV 




0.. 


.0.5 


7.44 


±0.75 


±0.62 


+1.04 
-0.72 


0.87 ±0.03 


1.01 ±0.02 


0.88 ±0.03 


0.5 


. . . 1 


12.1 


±1.0 


±0.4 


+1.3 
-1.2 


0.96 ±0.02 


1.03 ±0.02 


0.99 ±0.03 


1 . . 


. 1.5 


10.2 


±0.9 


±0.2 


+1.0 
-0.9 


0.99 ±0.02 


1.04 ±0.02 


1.02 ±0.03 


1.5 


. . . 2 


6.85 


±0.77 


±0.19 


+0.77 
-0.71 


1.00 ±0.02 


1.07 ±0.02 


1.08 ±0.02 


2 . . 


.2.5 


1.97 


±0.63 


±0.13 


+0.62 
-()..')() 


1.02 ±0.03 


1.07 ±0.04 


1.09 ±0.03 








52 < ^ 


f * < 75 GcV 


; 212 < W^.p 


< 285 GcV 




0.5 


. . .1 


1.17 


±0.33 


±0.08 


+0.07 
-0.11 


0.94 ±0.05 


1.02 ±0.05 


0.96 ±0.05 


1.. 


.1.5 


1.49 


±0.35 


±0.08 


+0.17 
-0.18 


0.99 ±0.04 


0.96 ±0.06 


0.96 ±0.07 


1.5 


. . .2 


1.24 


±0.32 


±0.12 


+0.19 
-0.15 


1.03 ±0.06 


1.02 ±0.05 


1.05 ±0.07 


2. . 


.2.5 


0.423 


±0.191 


±0.045 


+0.073 
-0.058 


1.02 ±0.08 


1.15 ±0.18 


1.17 ±0.18 



Table 4: Measured differential e^p cross section da/drj^^^ for inclusive jet photoproduction 
(Q^ < 1 GeV^j. The phase space of the measurement is divided into two regions ofWjp and 
three regions of E^^^ (see Tab. 1 caption for further details). 
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jyjet j-ange 


da/drj^''^ /^stat. ^syst. K.S. 


(1 + 5 frag) (1 + K.e) (1 + ^hadr) 




[nb] 





5 < E^^' < 12 GeV 



-1... 


-0.5 


16.2 


±0.4 


±1.1 


+2.0 
-1.8 


0.79 ±0.06 


1.29 ±0.06 


1.02 ±0.03 


-0.5 


...0 


18.1 


±0.4 


±1.4 


+2.1 
-2.0 


0.80 ±0.06 


1.45 ±0.08 


1.15 ±0.02 


0... 


0.5 


20.6 


±0.4 


±1.4 


+2.5 
-2. ,3 


0.73 ±0.06 


1.64 ±0.13 


1.19 ±0.02 


0.5. 


. .1 


24.1 


±0.5 


±1.7 


+2.7 
-2.8 


0.68 ±0.06 


1.89 ±0.27 


1.27 ±0.06 


1... 


1.5 


25.8 


±0.4 


±1.7 


+3.1 
-2.8 


0.69 ±0.06 


1.99 ±0.32 


1.35 ±0.10 



12 < Ei,"^ < 21 GeV 



-0.5 


...0 


0.584 


±0.081 


±0.034 


+0.063 
-0.102 


0.77 ±0.11 


1.12 ±0.07 


0.87 ±0.11 


0... 


0.5 


0.987 


±0.110 


±0.080 


+0.106 
-0.102 


0.91 ±0.09 


1.22 ±0.06 


1.11 ±0.11 


0.5. 


..1 


1.02 


±0.11 


±0.07 


+0.09 
-0.20 


0.89 ±0.09 


1.33 ±0.08 


1.19 ±0.09 


1... 


1.5 


0.803 


±0.096 


±0.060 


+0.175 
-0.110 


0.91 ±0.07 


1.53 ±0.14 


1.39 ±0.08 


1.5. 


..2 


1.12 


±0.11 


±0.08 


+0.26 
-0.15 


0.83 ±0.14 


1.74 ±0.33 


1.41 ±0.08 


2. . . 


2.5 


0.824 


±0.078 


±0.056 


+0.264 
-0.171 


0.88 ±0.09 


1.75 ±0.27 


1.53 ±0.14 



Table 5: Measured differential e+p cross section da/drj^^* for inclusive jet photoproduction 
(Q^ < 0.01 GcV^j, integrated over two E^'^^ ranges in the kinematic region 164 < W^p < 
242 GeV (see Tab. 1 caption for further details). 
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Ef' range {E^^') 


da /dE^^ Astat. ^syst. Ae.s. 


[GeV] 


[pb/GeV] 



164 < W^p < 242 GeV ; < q.oi GeV^ 



5. 


..8 


6.1 


39900 


±200 


±2600 


+4700 
-4400 


8.. 


.12 


9.3 


3840 


±60 


±250 


+690 
-570 


12. 


. .16 


14.3 


595 


±28 


±35 


+91 
-82 


16. 


. .21 


18.1 


126 


±12 


±7 


+33 
-18 



164 < W^p < 242 GeV ; < 1 GeV^ 



21. 


.28 


23.6 


31.3 


±0.4 


±0.8 


+3.0 
-2.9 


28. 


.35 


30.8 


6.66 


±0.19 


±0.18 


+0.56 
-0.57 


35. 


.42 


37.8 


1.73 


±0.10 


±0.06 


+0.19 
-0.16 


42. 


.52 


45.8 


0.415 


±0.042 


±0.014 


+0.047 
-0.044 


52. 


.62 


55.8 


0.0794 


±0.0179 


±0.0048 


+0.0131 
-0.0100 


62. 


.75 


66.4 


0.0143 


±0.0065 


±0.0004 


+0.0023 
-0.0021 



Table 6: Measured differential e~^p cross section da/dE^'^ for inclusive jet photoproduction 
integrated over —1 < t]^'^^ < 2.5 in the kinematic region 164 < Wyp < 242 GcV. Jets are 
defined using the cone algorithm with R = 1. For each range of E^^\ the average value of 
E^^* calculated with PYTHIA is given in the second column. The statistical ( Agtat.), systematic 
( Agyst.) and LAr energy scale ( Ag sJ uncertainties are shown separately. 



XT 


S{xt) 




Agyst. 


Ae... 


0.09 


0.0441 


±0.0011 


±0.0024 


±0.0078 


0.13 


0.0201 


±0.0014 


±0.0012 


±0.0028 


0.18 


0.00658 


±0.00106 


±0.00068 


±0.00241 


0.25 


0.00236 


±0.00005 


±0.00015 


±0.00022 


0.40 


0.000684 


±0.000053 


±0.000016 


±0.000073 


0.56 


0.000185 


±0.000047 


±0.000009 


±0.000028 



Table 7: Scaled 7p cross section at W^p = 200 GeV as a function of xt for \rf\ < 0.5. Jets 
are found with the cone algorithm (R = 1). The statistical (Agtat), systematic (Agyst.) and LAr 
energy scale fAg.sJ uncertainties are shown separately. 
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Figure 2: Top: differential e^p cross section for inclusive jet production as a function ofE^^^ in- 
tegrated over —1 < 7]^^^ < 2.5. The data are compared with LO and NLO QCD calculations 
using GRV photon PDFs and CTEQ5M proton PDFs. Bottom: relative difference between the 
data or different calculations and the NLO prediction with hadronisation corrections. The un- 
certainty associated with the LAr energy scale is shown as a hatched band. The shaded band 
displays the uncertainty on the NLO QCD predictions. The inside part shows the uncertainty 
associated with the hadronisation corrections, the outside part shows the uncertainty associ- 
ated with the choice of the renormalisation and factorisation scales and both uncertainties are 
added linearly. 
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Figure 3: Top: differential e^p cross section for inclusive jet production as a function ofE^^^ in- 
tegrated over -1 < r^-'^* < 2.5 /or 95 < W^p < 212 GeV (a) and 212 < W^p < 285 GeV 
(b). Bottom: relative difference between the data or different calculations and the NLO calcu- 
lation, including hadronisation corrections, based on GRV and CTEQ5M ( see Fig. 2 caption 
for further details). 
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Figure 4: Top: differential e^p cross section for inclusive jet production as a function ofE^^ in- 
tegrated over —1 < Tj^^^ < 2.5 and < 1 GeV^. The "low" E^^^ part, measured for 
< 0.01 GeV^, is corrected by a factor Rp which is the ratio of the photon fluxes in the 
two regions (see text). The photon-proton centre-of-mass energy is restricted to the range 



164 < W^p < 242 GeV. The data stemming from the analysis at "low" and "high" Ej, are 
indicated by empty and full points respectively. Bottom: relative difference between the data or 
LO QCD prediction and the NLO calculation, including hadronisation corrections, based on 
GRV and CTEQ5M (see Fig. 2 caption for further details). 
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Figure 5: Top: differential e^p cross section for inclusive jet production as a function ofrj^^^ in- 
tegrated over 21 < E^^^ < 75 GeV. Bottom: relative difference between the data or different 
calculations and the NLO calculation, including hadronisation corrections, based on GRV and 
CTEQ5M (see Fig. 2 caption for further details). 
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by using GRV photon PDFs and CTEQ5M proton PDFs (see Fig. 2 caption for further details). 
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Figure 7: Differential e^p cross section for inclusive jet production as a function ofrj^^^ inte- 
grated over various E^^^ and W^p ranges. The data are compared with NLO QCD predictions 
obtained by using different photon and proton PDFs (see Fig. 2 caption for further details). 
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Figure 8: Differential e^p cross section for inclusive jet production as a function of 7]^^^ inte- 
grated over various E^^^ ranges. The data are compared with LO and NLO QCD predictions 
obtained by using GRV or AFG photon PDFs and CTEQ5M proton PDFs (see Fig. 2 caption 
for further details). 
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Figure 9: Scaled 'jp cross section at W^p = 200 GeV for inclusive jet production as a function 
ofxxfor |?7*| < 0.5. Jets are found with the cone algorithm (R = 1). The data are compared 
with measurements from UAl [4, 5] and DO [6, 7] of inclusive jet production in pp collisions at 
various cms energies. The predictions of PYTHIA for •yp and for the resolved photon contribu- 
tion are also shown, multiplied by a factor 1.2 to match the normalisation of the data. 
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